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LIST OF ABBREVIATIONS 



ky - surface area of water droplets per volume of cooling tower 

C p> - constant pressure specific heat of moist air 

C pw - constant pressure specific heat of liquid water 

C, - derivative of saturation air enthalpy with respect to 

temperature 

h 4 - enthalpy of moist air per mass of dry air 

h c - convection heat transfer coefficient 

h D - diffusion mass transfer coefficient 

h„ - enthalpy of water above reference state for liquid water 

at T rcf 

h s - enthalpy of saturated air 

xiV - mass flow rate of water 

m a - mass flow rate of air 

m* - ratio of air flow rate to water effective capacitance rate 

for effectiveness analysis 

Le - Lewis number 

Ntu - overall number of transfer units 

Q - actual heat transfer rate 

T ref - reference temperature for zero enthalpy of liquid water 
(32 °F) 

T w - water temperature 

V - volume 

e, - air-side heat transfer effectiveness 

w, - air humidity ratio 

w s - humidity ratio of saturated air 
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Additional Subscripts 



a - air stream conditions 

i - inlet state conditions 

o - outlet state conditions 

T - total 

w - water stream conditions 
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SUMMARY 



Optimization of the design of chiller/cooling tower systems 
for large office buildings require an accurate detailed model for 
the cooling tower. This model must be precise and allow for 
varying air and water flow rates, ambient air conditions and 
entering water temperatures. The optimization computational 
requirement is very time consuming, however, due to the tower model 
requiring solution of three simultaneous differential equations. 

The purpose of this project was to develop a computer program 
utilizing an effectiveness model to reduce computational time while 
maintaining good correlation to existing programs. The methodology 
involved the use of a tower model program developed by Eric Weber 
[1988] and the modifications made unto by Thane Joyce [1990] that 
integrated a centrifugal chiller model, cooling tower model, and 
system optimization strategy. The effectiveness model generated by 
James Braun [1988] was utilized in place of Weber's tower model. 

The results are very encouraging. The computational time 
required to optimize the design was reduce approximately 75% while 
maintaining an accuracy of tower water temperature range of 
approximately 15%. This reduction in computer time will allow an 
engineer the flexibility to examine many systems and determine the 
optimum design in a shorter time frame. 
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INTRODUCTION 



The design of heating, ventilation and air conditioning (HVAC) 
for large buildings tasks the engineer with manipulating a plethora 
of design parameters. With limitations on funding and typically a 
short fused deadline, engineers rely on rule-of-thumb or past 
experience when designing these systems. No attempt is made to 
optimize the system as a unit bringing together all component 
parts: pumps, chiller, cooling tower, pipe size, etc. 

An optimization methodology for the chiller/cooling tower 
system was developed by Thane Joyce. The results of his work show 
that an optimized system can save both construction and operating 
costs. The computational requirements to perform the optimization, 
however, is very time consuming and therefore the methodology was 
not applicable as a design tool. In order to make the program more 
useful, the computational time needed to be reduced. This is 
accomplished by employing Braun's effectiveness model in lieu of 
the Weber cooling model in Joyce's program. The resulting program 
was then compared to the original Joyce program for accuracy and 
reduction in computational time. 
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DEVELOPMENT OF EFFECTIVENESS COMPUTER MODEL 



The effectiveness model used to replace the tower model 
developed by Weber is based on the work by James Braun. A 
schematic of a counterflow cooling tower is provided in figure 1. 
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Figure 1. Schematic of a Counterflow Cooling Tower [Braun, 
1988 ] 



2 



The following assumptions were made by Braun to derive the basic 
model equations: 

1. Heat and mass transfer in the direction normal to flows only. 

2. Negligible heat and mass transfer through tower walls to 
environment. 

3. Negligible heat transfer from the tower fans to the air or 
water streams. 

4. Constant water and dry air specific heats. 

5. The mass fraction of water vapor in the mixture of air and 
vapor is approximately equal to the humidity ratio. 

6. Uniform temperature throughout the water stream at any cross 
section. 

7. Uniform cross-sectional area of the tower. 

From steady-state energy and mass balances on an incremental 
volume, refer to figure 1, the following differential equations are 
developed by Braun and Weber: 
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where 





Knowing the inlet conditions, the number of transfer units, and the 
Lewis number, the exit state of both the air and water streams can 
be determined. Using these same equations, Weber took an iterative 
approach with respect to the outlet humidity ratio and the exit 
water temperature, numerically integrating the equations over the 
entire volume from air inlet to outlet. 

The effectiveness model simplifies this analysis considerably. 
By employing Merkel's assumptions of neglecting water loss and 
setting the Lewis number equal to unity, the three differential 
equations reduce to: 



Equation (5) can be written in terms of air enthalpies only by 
introducing a new variable, C,, termed the saturation specific heat: 
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where 



<?, = [ 



dh, 

dT 
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r.r. 



(7) 



C, is equivalent to the derivative of the saturated air enthalpy 
with respect to temperature evaluated at the water temperature and 
has the units of specific heat. Braun indicates that if the 
saturation air enthalpy were linear with respect to temperature 
then the exit conditions could be solved for analytically. His 
figure 2, shown below, depicts the variation of the saturation 
enthalpy with temperature for typical water inlet and outlet points 
along with a straight line connecting the points. 




Figure 2. Saturation Air Enthalpy versus Temperature [Braun, 
1988] 
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It is clear form the figure that the saturation air enthalpy does 
not vary linearly with temperature, however, by choosing an 
appropriate average slope between the inlet and outlet states, an 
air-side effectiveness relationship in terms of C, can be developed. 
This air-side effectiveness is defined as the ratio of the actual 
heat transfer to the maximum possible heat transfer if the exiting 
air were saturated at the temperature of the incoming water (i.e., 

Writing the heat transfer in terms of this effectiveness 
yields: 

Q = e 0 m a (h swi - hj ( 8 ) 



where 



e a = air-side heat transfer effectiveness 
1 - exp ( - Ntu (1 - m * ) ) 

1 - m * exp ( - Ntu (l - m * ) ) 



( 9 ) 



and 



i(Cpw/C s ) 



( 10 ) 



The average value for the saturation specific heat is estimated as 
the average slope between the inlet and outlet state points: 
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The outlet air enthalpy and water temperature can then be 
calculated using the following equations: 
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In the computer model employing the air-side effectiveness, the 
water loss is neglected, so m^; is set equal to m^. Since C s 
depends upon T wo , the solution for the exit conditions is iterative 
using equations 9 through 13. The saturation specific heat, C,, is 
only weakly dependent upon T wo , however, so any reasonable initial 
guess for T wo (such as the inlet wet bulb temperature of the air) 
typically results in convergence in only two iterations. 
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CONCLUSIONS 



The results of this exercise are very promising. The 
computational time savings for the effectiveness method was as much 
as 83% over the Joyce method, while the percent error in 
calculating the energy consumed was less than 14% for all cases. 

The first portion of the project is dedicated to evaluating 
the affect of employing the effectiveness model on the delta 
temperature across a counterflow cooling tower. Weber's computer 
program is used to generate the control data set. The revised 
program generated the test data set. Figure 3 displays the results 
o f this evaluation. Numerical data is found in Appendix B. 

ACCURACY OF WEBER vs EFFECTIVENESS 
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Figure 3. Accuracy of Weber versus Effectiveness Model 
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The effectiveness model computer program is provided as Appendix A. 
This model gives a larger delta temperature across the cooling 
tower in all cases compared with the exact data set. One reason 
why the effectiveness model overpredicts the tower delta 
temperature is that Weber's model is developed for a crossflow 
cooling tower while Braun's analysis is for a counterflow cooling 
tower. If improved accuracy in calculating the delta temperature 
across the tower is desired, a crossflow correction factor can be 
applied to the effectiveness model. Linearization of saturation 
enthalpy between the inlet and outlet points tends to overpredict 
the tower heat transfer and is another factor contributing to this 
error . 

To test and compare the effectiveness model to the Weber 
cooling tower model in Joyce's program, several variables are 
manipulated while holding others constant. The variables 
investigated are: 1) air flow rate through tower (cfm) , 2) tower 
size, and 3) pump size (gpm) . These variables are considered the 
primary parameters for designing cooling towers. Figures 4 through 
9 show the results of the comparisons for these various parameters. 
The results are encouraging and in each case the trends for both 
models are similar. Appendix C contains the numerical data used in 
generating these figures. 

Figure 4 depicts the plot of effectiveness data versus Joyce 
data for various air flow rates through the cooling tower (CFM) . 
The trends in both computational methods for this design parameter 
are very similar. 
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ACCURACY OF JOYCE VS EFFECTIVENESS 

FOR VARIOUS CFM'S 
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Figure 4 . Plot of Effectiveness versus Joyce Data for Various 
CFMs. 



The time savings for this variable are substantial, as shown in 
figure 5. 

Six different tower air flow rates were used while holding the 
other variables constant. The maximum error, of 12.8%, occurs at 
a flow rate of 160000 CFM. The time savings are between 33 and 37 
minutes, or roughly 81% improvement over the more exact Joyce 
model . 

The next variable manipulated is tower size. Actual 
manufacture's performance data for Baltimore Aircoil series 3000 is 
used as input parameters. The results are very favorable giving 
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Figure 5. Time Savings at Various CFMs. 



similar trends for both models, figure 6. The time savings using 
the effectiveness model over the Joyce method for this parameter 
are also quite impressive, refer to figure 7. 

Six different tower sizes were used in the analysis. The 
effectiveness model underpredicts the tower fan energy used in all 
cases. The largest error, 13.7%, is found at tower number 19 which 
denotes a Baltimore Aircoil, series 3000, model number 3269. 
Savings in computational time is very impressive for this test 
variable. Approximately 34 minutes of saved time is seen for all 
cases, roughly an 80% improvement in computational run time. 
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ACCURACY OF JOYCE VS EFFECTIVENESS 

FOR VARIOUS TOWER SIZES 




Figure 6. Plot of Effectiveness versus Joyce Data for Various 
Tower Sizes. 



The final variable manipulated is condenser water pump size. Pump 
sizes between 750 GPM and 1500 GPM were used and the affect on both 
fan and pump energy (KWH) were tabulated. The results are 
promising, producing like trends for both the effectiveness model 
and Joyce method, refer to figure 8. Computational run time 
savings are also favorable for this parameter. Figure 9 presents 
the time savings for the new program as compared with the existing 
program for various pump sizes. 
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Time Savings at Various Tower Sizes 



Again the effectiveness model underpredicts the fan and pump energy 
as compared to the Joyce method in all cases. The trends for both 
models, however, are very similar giving the same minimum and 
maximum point. The maximum error is roughly 10% and occurs at a 
pump size of 1000 GPM. Computer run time savings is very 
impressive at approximately 34 minutes or 81% as compared to the 
control data program. 

As is shown in all the figures, the trends for both the 
effectiveness model and the Weber model in Joyce's program are very 
similar. 
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ACCURACY OF JOYCE VS EFFECTIVENESS 

PUMP SIZE 




Figure 8. Plot of Effectiveness versus Joyce Data for Various 
Pump Sizes. 



In order to improve the accuracy of the effectiveness method, 
crossflow correction factor should be applied and the water loss 
from the cooling tower should be included. The affect on computer 
run time, however, will be adverse as several other iteration loops 
would have to be utilized. Another improvement could be using a 
more accurate curve fit for the saturation specific heat C t . This 
will also take additional computational time and would result in 
only a minimal improvement in accuracy. 
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Figure 9. Time Savings at Various Pump Sizes, 



The results of this project are very encouraging. For a 
general trend and a first time run through finding the optimum 
design this revised effectiveness model fits the bill. The 
reduction in computational time will allow an engineer flexibility 
in investigating various systems reaching an optimum design in a 
shorter time frame. 
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APPENDIX A 



EFFECTIVENESS COMPUTER PROGRAM 
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/ 



EFFECTIVENESS PROGRAM MODEL 



' 15 OCT 91 

' THIS PROGRAM CALCULATES THE ANNUAL ENERGY COST, AND THE 
' CAPITAL CONSTRUCTION COST OF THE CONDENSER SIDE OF A CENTRIFUGAL 
' CHILLER PLANT. VARIES THE CONDENSER GPM TO MAINTAIN OPT GPM FOR 
' A PARTICULAR TEMP. BIN ASSUMING 1071 GPM MAX. USED TO COMPARE 
' FIXED SPEED DESIGN WITH VARIABLE SPEED. 

/ 

' INPUT PARAMETERS : 

t 

' 1) THE CITY (ATLANTA=1, L.A.=2, N.Y.=3) 

' 2) THE COOLING TOWER NO. CORRESPONDING TO ONE OF 

' THE COOLING TOWERS IN THE DOCUMENTATION. 

' 3) THE CONDENSER WATER PIPE SIZE IN INCHES. 

' 4) THE COOLING TOWER CFM. 

' 5) THE NO. OF CONDENSER PASSES IN THE CHILLER. 

' 6) THE MAXIMUM FLOWRATE (GPM) FOR VARIABLE SPEED 

' PUMP. 



' **** THIS PROGRAM IS A MODIFICATION OF THE PROGRAM ANUALFIX 

' WRITTEN BY ERIC WEBER IN CONJUNCTION WITH HIS MASTERS 

' THESIS ENTITLED "MODELING AND GENERALIZED OPTIMIZATION OF 

' COMMERCIAL BUILDING CHILLER/ COOLING TOWER SYSTEMS" 

I 

' UPDATES 

' LAST UPDATE TO ANUALFIX BY EDW ON 23 OCT 1988 

' MOD TO PREVENT RECALC OF PUMP HP BY CTJ ON 3 SEPT 89 

' MOD TO PREVENT RECALC OF FAN HP BY CTJ ON 5 SEPT 89 

' MOD TO TWR CALC. FOR COOLING WITH FAN OFF BY CTJ ON 15 SEP 89 

' MODI TO BISECTION USE OLD CW TEMP FOR BRACKETS BY CTJ ON 29 

' SEP 89 

' MOD 2 TO BISECTION TO CHECK SIGNS ON BRACKETS BY CTJ ON 30 

' SEP 89 

' MOD 3 TO BISECTION WARN IF SIGNS ARE SAME BY CTJ ON 1 OCT 89 

' MOD: STD WTHR DATA NOT COMPR. WTHR DATA BY CTJ ON 15 OCT 89 

' MOD TO ADD WTHR DATA FOR N.Y. AND L. A. BY CTJ ON 20 OCT 89 

' MOD TO ADD ENERGY BILL CALC. FOR 3 CITIES BY CTJ ON 23 OCT 89 

' MOD TO ADD CAPITAL COST FOR 3 CITIES BY CTJ ON 25 NOCT 89 

' MOD TO USE BRENT METHOD VICE OF BISECTION BY CTJ ON 27 OCT 89 
' MOD TO ASSIGN TCI 70 IF PREVIOUS HIGHER TDB ASSIGNED TCI TO 

' 70, AND CHECK FOR INCREASE IN DELTATCOND BY CTJ ON 29 OCT 89 

' MOD TO ADD VARIABLE SPEED PUMP BY CTJ ON 29 OCT 89 

' UPDATE TO INCLUDE THE EFFECTIVENESS MODEL TO CALCULATE THE 

' OUTPUT TEMP OF THE COOLING TOWER ON 15 OCT 91 BY JDR 

/ 

DECLARE SUB EFFECTIVE (NTUi , TWIN!, MDOTW! , MDOTA! , TWB! , TDB!, 
EXPTWO ! , TWOUT ! ) 
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DECLARE SUB BRENT (TWIN1!, TWIN2 I , DIFF1!, DIFF2 ! , TW0UT1 ! , 
TWOUT2 ! , NTU! , f da ! , MDOTW! , HDOTA! , TDB! , TWB! , DELTATCONDl! , 
CTR2%, TCI!, TC2 ! , CITY!) 

DECLARE SUB SYSCOST (PMPHP! , FANHP! , TWRCOST!, diameter!, CITY!, 
CAPCST ! ) 

DECLARE SUB ENRGYCST (KWPK!(), KWHTOT! () , BILLDEM! ( ) , BILL!(), 
YRBILL! , PMPKWH ! ( ) , FANKWH! () , CMPKWH!(), CITY!) 

DECLARE SUB Hairm (T2!, Tl!, W2 ! , OMEGMIXTURE! , HAIRMIX!) 

DECLARE SUB Hairs (TF! , WS!, HSAT ! ) 

DECLARE SUB FLOW (GPM! , CFM! , TDB!, MDOTW!, MDOTA ! ) 

DIM HRSDAT(20) , KWPK(20), KWHTOT(20), BILLDEM (20) , BILL(20), 

TWBR(20) , TDBR(20) , PMPKWH(20), FANKWH ( 20 ) , CMPKWH(20) 

I 

' *************************** 

' USER INPUT 

/ 

INPUT "INPUT CITY NUMBER (1-3) ", CITY 

INPUT "INPUT TOWER NUMBER (1-19) ", i 

INPUT "INPUT DESIRED CONDENSER WATER PIPE SIZE IN INCHES ", 
diameter 

INPUT "INPUT DESIRED TOWER CFM ", CFM 

INPUT "INPUT NUMBER OF CONDENSER PASSES ", PASSES 

INPUT "INPUT THE MAXIMUM GPM FOR A VARIABLE SPEED PUMP ", MAXGPM 
****************************************** 

'STORE THE START TIME FOR THE ITERATIONS 
'***************************************** 

/ 

STARTIME$ = TIME$ 

********************************************************** 

' GET THE TOWER CONSTANT AND THE EXPONENT FROM DATA FILE 
• ********************************************************* 



KWHSUM = 0! 

TWR = i 

OPEN "b: tower.dat" FOR INPUT AS #1 
DO WHILE i > 0 

INPUT #1, TNUM, EXPONENT, CONSTANT, GPML, GPMH, XXX, FHPCA , 
FHPCB, TWRCST 
i = i - 1! 

LOOP 

CLOSE #1 
/ 

' THIS SECTION SETS UP THE DRY BULB TEMPERATURE BINS 
'******************************************************* 
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FOR k = 1 TO 10 

TDBR(k) = 102 - ( (k - 1) * 5) 

NEXT k 
/ 

***************************************************************** 
' THIS SECTION GETS THE WEATHER DATA FOR THE DRY BULB TEMP. BINS ' 
' FOR EACH MONTH — JAN(i=l) thru Dec(i=12) 
***************************************************************** 
/ 

FOR i = 1 TO 12 
/ 

'********************************** 

' INITIALIZE MONTHLY KWH AND PEAKKW 
*********************************** 

KWHTOT(i) = 0! 

KWPK(i) = 0! 

PMPKWH(i) = 0! 

FANKWH(i) = 0! 

CMPKWH(i) = 0! 

CWSTEMPCHK = 0 

CTR% = 0 

HRSMNTH = 0! 

SELECT CASE CITY 
CASE 1 

OPEN "b: atlwth.dat" FOR INPUT AS #3 
CASE 2 

OPEN "b: lawth.dat" FOR INPUT AS #3 
CASE 3 

OPEN "b: nywth.dat" FOR INPUT AS #3 
END SELECT 
k = i 

DO WHILE k > 0! 

FOR m = 1 TO 10 

INPUT #3, HRSDAT(m), HOURSSUM, TWBR(m) 

NEXT m 
k = k - 1 
LOOP 

CLOSE #3 

/ 

9 ****************************** 

'SELECT THE WEATHER DATA BIN 
9 ****************************** 

/ 

FOR j = 1 TO 10 
TWB = TWBR(j) 

EXPTWO = TWB 



'WET BULB TEMP 

'EXPECTED TEMP WATER FROM TOWER 
' FOR EFF MODEL 



'MONTHLY KWH 

'MAX MONTHLY KW FOR AT LEAST 15% OF HRS 
' IN MONTH 



'CHECK VALUE TO INSURE CWS TEMP DOESNT 
' RISE WITH DECREASING TDB, AND LOAD 
'TRACKS # OF TEMP BINS WITH NON-O TWB 
' PER MONTH 
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TDB = TDBR(j) 


'DRY BULB TEMP 




HOURS = HRSDAT(j) 


'MONTHLY HOURS AT TDB & TWB 
' 6AM TO 6PM 


FROM 


IF TWB = 0 THEN GOTO 
CTR% = CTR% + 1 


GETNEXTWEATHER 




CTR1% = 0 


'TRACKS # OF ITTERATIONS 
' RESTARTA 


THRU 


CTR2% = 0 


'TRACKS # OF ITTERATIONS 
' BRENT MTHD 


THRU 


CTR3% = 0 


'TRACKS # OF ITTERATIONS 
' TRYAGAIN 


THRU 


CTR4% = 0 


'TRACKS # OF ITTERATIONS 
' CYCLEFAN 


THRU 



'************************************** 

' CALCULATION OF THE EVAPORATING TEMP. 
'************************************** 
/ 

EVAPGPM = 1200 
TE2 = 44 



'****************************** 

' EVAP LOAD ALGORITHM 
• ****************************** 

/ 

GPM = -1242 + (93.90529 * TDB) - (1.45875 * TDB A 2) + (.007462 

* TDB A 3) 

IF GPM <= 775 THEN GPM =750 
tons = -531.578 + 10.526 * TDB 
IF tons >500 THEN 
tons = 500 
END IF 

IF tons < 100 THEN 
tons = 100 
END IF 

QDOTEVAP = tons * 12000 

DELTATEVAP = QDOTEVAP / (EVAPGPM * 499.27) 

UAEVAP = 1000 * (134.11603# + 1284.3418# * (tons / 500) - 

640.54974# * ((tons / 500) * 2)) 

LMTDevap = QDOTEVAP / UAEVAP 

TE = (-TE2 - DELTATEVAP + TE2 * EXP (DELTATEVAP / LMTDevap)) / 
(EXP (DELTATEVAP / LMTDevap) - l!) 

t 

f *** ************** ********** ******** 

' ESTIMATE OF THE KW FOR COMPRESSOR 

• *********************************** 

! 

KW = 30 + .61 * (tons) 



' FLOWRATE THRU EVAPORATOR 
'FIXED CHILLED WATER SUPPLY 
' TEMP. 
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TWIN1 = 50.5 
TWIN2 = 110 

/ 

'***************************************************** 

' CALCULATE THE HEAT BURDEN FOR THE CHILLER CONDENSER 
'***************************************************** 

/ 

RESTARTA: 

QDOTCOND = QDOTEVAP + (KW * 3413) 

DELTATCOND1 = QDOTCOND / (GPM * 499.27) 

CALL FLOW (GPM, CFM, TDB, MDOTW, MDOTA) 

fda = MDOTW * CONSTANT * ((MDOTW / MDOTA) A EXPONENT) 

NTU = CONSTANT * ((MDOTW / MDOTA) A (EXPONENT +1)) 

/ 

' ******************************* 

' BRACKETING FOR BRENT METHOD 
'******************************** 

/ 

CALL EFFECTIVE (NTU, TWIN1, MDOTW, MDOTA, TWB, TDB, EXPTWO, 

TWOUT1) 

DELTA1 = TWIN1 - TWOUT1 
DIFF1 = DELTATCOND1 - DELTA 1 

CALL EFFECTIVE (NTU, TWIN2 , MDOTW, MDOTA, TWB, TDB, EXPTWO, 

TWOUT2 ) 

DELTA2 = TWIN2 - TWOUT2 
DIFF2 = DELTATCOND1 - DELTA2 
IF SGN(DIFFl) = SGN (DIFF2 ) THEN 
TWIN1 = TWIN1 - 5 
TWIN2 = TWIN2 + 5 
SELECT CASE CITY 
CASE 1 

OPEN "B: ATLEFF.DAT" FOR APPEND AS #4 
CASE 2 

OPEN "B: LAEFF.DAT" FOR APPEND AS #4 
CASE 3 

OPEN "B: NYEFF.DAT" FOR APPEND AS #4 
END SELECT 

PRINT #4, "ROOT NOT BRACKETED FOR BRENT. PROG ASSIGNING NEW 
BRACKETS" 

PRINT #4, "MONTH= i, "BIN =" ; TDB 
CLOSE #4 

PRINT "ROOT NOT BRACKETED FOR BRENT. PROG ASSIGNING NEW 
BRACKETS" 

PRINT "MONTH= "; i, "BIN =" ; TDB 
GOTO RESTARTA: 

END IF 

CTR1% = CTR1% + 1 

CALL BRENT (TWIN1, TWIN2 , DIFF1 , DIFF2 , TWOUT1, TWOUT2 , NTU, fda, 
MDOTW, MDOTA, TDB, TWB, DELTATCOND1 , CTR2%, TCI, TC2 , CITY) 

DELTA 3 = TWIN2 - TWOUT2 
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/ 



************************************************** 

' DETERMINATION OF CONDENSER PERFORMANCE (UA,Tc) 

9 ************************************************* 

/ 

GPMCOND = GPM 

UAcond = 608.5483400000001# + 130.64856# * (tons / 500) - 
192.20566# * ((tons / 500) A 2) 

UAcond = UAcond + 1203.80505# * (GPM / (PASSES * 1500)) - 
557.4549# * ((GPM / (PASSES * 1500)) A 2) 

UAcond = UAcond * 1000 
LMTDCOND = QDOTCOND / UAcond 

TC = (TCI - (TCI + (TC2 - TCI)) * EXP((TC2 - TCI) / LMTDCOND)) / 
(1! - EXP ( (TC2 - TCI) / LMTDCOND)) 

'*************************************************************** 

'DETERMINE MODEL KW USING THE CARNOT KW/TON AND TWO EFFICIENCIES 
'*************************************************************** 

/ 

KWCARNOT = (((TC + 460) / (TE + 460)) - II) * 3.517 
NfTON = 4.5869 * (tons / 500) - 8.16536 * ((tons / 500) A 2) + 
6.65014 * ((tons / 500) A 3) - 2.07167 * ((tons / 500) A 4) 

NfISENT = .516061 + .31889 * ((KWCARNOT - .3) / .2) - .14911# * 
(((KWCARNOT - .3) / .2) A 2) 

KWTONMODEL = KWCARNOT / (NfISENT * NfTON) 

KWTONGUESS = KW / tons 

/ 

^ ********************************************** 

' COMPARE MODEL KW/TON TO THE GUESSED KW/TON 
9 ********************************************** 

IF ABS (KWTONMODEL - KWTONGUESS) > .0001 THEN 
KW = KWTONMODEL * tons 
GOTO RESTARTA 
END IF 
TCI A = TCI 
TC2A = TC2 

/ 

'**************************************************************** 
' IF THE ENTERING COND. WATER TEMP. IS BELOW 80 THEN FIX IT AT 80 
' AND CYCLE THE COOLING TOWER FAN. 
'**************************************************************** 



IF TCI < 80 THEN 
GOTO CYCLEFAN 
END IF 

IF CTR% > 1 THEN 
GOTO CALCTOTALKW 
END IF 

FANHP = EXP ( FHPCA + ( (LOG (CFM) - 10.5) / 21) * FHPCB) 
FanKW = ((FANHP * .7457) / .85) 

GOTO CALCTOTALKW 
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/ 

'**************************************************************** 



' THIS SECTION WILL HOLD CONDENSER WATER INPUT TEMPERATURE TO 
' 80 F AND CYCLE THE TOWER FAN USING THE TOWER TEMPERATURE DROP 
' VERSUS THE CONDENSER TEMPERATURE DROP 

r 

> **************************************************************** 



CYCLEFAN: 

TRY AGAIN: 

CTR3% = CTR3% + 1 

QDOTCOND = QDOTEVAP + (KW * 3413) 

DELTATCOND = QDOTCOND / (GPM * 499.27) 

CALL FLOW (GPM, CFM, TDB , MDOTW, MDOTA) 

fda = MDOTW * CONSTANT * ( (MDOTW / MDOTA) A EXPONENT) 

NTU = CONSTANT * ((MDOTW / MDOTA) A (EXPONENT +1)) 

TCI = 80 

TC2 = TCI + DELTATCOND 

CALL EFFECTIVE (NTU, TC2 , MDOTW, MDOTA, TWB, TDB, EXPTWO, TWOUT1) 
DELTATTOWER = TC2 - TWOUT1 

UAcond = 608.5483400000001# + 130.64856# * (tons / 500) - 
192.20566# * ((tons / 500) A 2) 

UAcond = UAcond + 1203.80505# * (GPM / (PASSES * 1500)) - 
557.4549# * ((GPM / (PASSES * 1500)) A 2) 

UAcond = UAcond * 1000 
LMTDCOND = QDOTCOND / UAcond 

TC = (TCI - (TCI + (TC2 - TCI)) * EXP((TC2 - TCI) / LMTDCOND)) / 
(1! - EXP ( (TC2 - TCI) / LMTDCOND)) 

KW CARNOT = ( ( (TC + 460) / (TE + 460)) - l!) * 3.517 
NfTON = 4.5869 * (tons / 500) - 8.16536 * ((tons / 500) * 2) + 

6.65014 * ((tons / 500) A 3) - 2.07167 * ((tons / 500) A 4) 

NfISENT = .516061 + .31889 * ( (KWCARNOT - .3) / .2) - .14911# * 
( ( (KWCARNOT - .3) / .2) A 2) 

KWTONMODEL = KWCARNOT / (NfISENT * NfTON) 

KWTONGUESS = KW / tons 

IF ABS (KWTONMODEL - KWTONGUESS) > .0001 THEN 
KW = KWTONMODEL * tons 
GOTO TRYAGAIN 
END IF 

FANHP = EXP (FHPCA + ((LOG (CFM) - 10.5) / 2!) * FHPCB) 

TOWERDUTY = (DELTATCOND - (.1 * DELTATTOWER)) / (DELTATTOWER * 
•9) 

FanKW = ((FANHP * .7457) / .85) * TOWERDUTY 
CTR4% = CTR4% + 1 

/ 

'**************************************************************** 

' CALCULATION OF THE POWER REQUIRED BY THE CONDENSER WATER PUMP 
, **************************************************************** 
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CALCTOTALKW: 

PipeD = diameter / 12 

Condenservelocity = (GPM * .002228) * (PASSES / 1.15241) 
Pipevelocity = (GPM * .002228) / ((3.1416 * (PipeD A 2) ) / 4) 
REd = (Pipevelocity * PipeD) / 8.64E-06 
IF REd < 2300 THEN 
pipefrict = 64 / REd 
ELSE 

• ********************************************* 

' FROM EQU. 6.64a IN FLUID MECHANICS BY WHITE 

9 ********************************************* 

AA = (6.9 / REd) + ((.00015 / (PipeD * 3.7)) A 1.11) 

BB = LOG (AA) / LOG (10!) 
pipefrict = (l! / (-1.8 * BB) ) A 2 
END IF 

PHtower = 11 
SELECT CASE PASSES 
CASE 1 

PHcond = EXP(-13. 95945 + 20.21629 * (LOG (GPM) / 10)) 

CASE 2 

PHcond = EXP(-10. 96759 + 19.02989 * (LOG (GPM) / 10)) 

CASE 3 

PHcond = EXP (-8. 96627 + 17.76877 * (LOG (GPM) / 10)) 

CASE 4 

PHcond = EXP (-8. 76329 + 18.62906 * (LOG (GPM) / 10!)) 

END SELECT 

PipeL = 200 

PipefL = 450 * PipeD 

PipestrainL = 250 * PipeD 

Ltotal = PipeL + PipefL + PipestrainL 

PHpipe = ((Pipevelocity A 2) / 64.4) * ((pipefrict * Ltotal) / 

PipeD) 

PHTOTAL = PHcond + PHtower + PHpipe 
PERGPM = (GPM / MAXGPM) * 100 

PMPEFF = 10.57993 + (2.436639 * PERGPM) - (.01966 * (PERGPM A 2)) 
VFDEFF = 4.983855 + (1.846339 * PERGPM) - (.01029 * (PERGPM A 2)) 
PMPHP = (GPM * PHTOTAL) / (3960 * PMPEFF / 100) 

PumpKW = (PMPHP * .7457) / (VFDEFF / 100) 

/ 

'******************************************************* 

' CALCULATION OF THE TOTAL POWER REQUIRED BY THE SYSTEM 
'******************************************************* 



TOTALKW = PumpKW + FanKW + KW 'TOTAL SYSTEM KW 
KWHBIN = TOTALKW * HOURS 'ENERGY CONSUMPTION OF BIN 

KWHTOT(i) = KWHTOT(i) + KWHBIN 'ENERGY CONSUMPTION OF MONTH 
KWHSUM = KWHSUM + KWHBIN 'ENERGY CONSUMPTION OVER YEAR 

PKHOURS = .15 * HOURS SUM 

IF (HOURS >= PKHOURS) AND (TOTALKW > KWPK(i)) THEN KWPK(i) = 
TOTALKW 
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HRSMNTH = HRSMNTH + HOURS 

PMPKWH(i) = PMPKWH(i) + (PumpKW * HOURS) 

FANKWH(i) = FANKWH(i) + (FanKW * HOURS) 

CMPKWH(i) = CMPKWH(i) + (KW * HOURS) 

/ 

************************************ 

'STORE THE END TIME FOR ITERATIONS 
/ *********************************** 
t 

ENDTIME$ = TIME$ 

/ 

'**************************************************************** 
7 PRINTING OF THE RESULTS 

***************************************************************** 



SELECT CASE CITY 
CASE 1 

OPEN "B: ATLEFF.DAT" FOR APPEND AS #4 
CASE 2 

OPEN "B: LAEFF.DAT" FOR APPEND AS #4 
CASE 3 

OPEN "B:NYEFF. DAT" FOR APPEND AS #4 
END SELECT 

PRINT #4, "TIME START = STARTIME$ ; "TIME END = ", ENDTIME$ 
PRINT #4, "ctrl=" ; CTR1%, "ctr2 = "; CTR2%, " ctr3 = CTR3%; "" 

PRINT #4, "I =" } i, " J=" ; j 

PRINT #4 \ "TDB (DEG F) ’ TWB (DEG F) HOURS 
PRINT #4, ; TDB, TWB, HOURS, TWR 

PRINT #4, "GPM %GPM CFM 

PRINT #4, ; GPM, PERGPM, CFM, tons 
IF TCI >80 THEN 

PRINT #4, "CWS (DEG F) CWR(DEG F) DELTATCOND 

PRINT #4, ; TC1A, TC2A, DELTATCOND 1 , DELTA3 

ELSE 

PRINT #4, "CWS (DEG F) CWR(DEG F) DELTATCOND 

TOWERDUTY" 

PRINT #4, ; TCI, TC2 , DELTATCOND, DELTATTOWER, TOWERDUTY 

END IF 

PRINT #4, "PIPE D (FT) Pipevelocity (FT/SEC) REd pipe 

frict " 



TOWER " 
TONS" 

DELTATTOWER" 

DELTATTOWER 



PRINT #4, ; PipeD, Pipevelocity, REd, pipefrict 
PRINT #4, "pipeL total PHpipe PHcond PHtower 

PHtotal(FT Water)" 

PRINT #4, ; Ltotal, PHpipe, PHcond, PHtower, PHTOTAL 

PRINT #4, "PASSES Condenservelocity" 

PRINT #4, ; PASSES, Condenservelocity 
/ 



PRINT #4, "UACOND LMTDCOND NfISENT 

COND.TEMP" 

PRINT #4, ; UAcond , LMTDCOND, NfISENT, NfTON, TC 

PRINT #4, "KWTONGUESS KWTONMODEL KWHBIN 



NFTON 



KWHTOT" 
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PRINT 


#4, 


; KWTONGUESS , 


KWTONMODEL, KWHBIN, KWHTOT (i) 


PRINT 


#4, 


"HOURS 


HOURSSUM" 


PRINT 


#4, 


; HOURS, HOURSSUM 


PRINT 


#4, 


" CHILL ERKW 


PUMPKW FANKW TOTALKW 


PRINT 


#4, 


; KW , PumpKW , 


FanKW, TOTALKW 


PRINT 


#4, 


• • It 




PRINT 


#4, 


19 It 




PRINT 


#4, 


91 99 




PRINT 


#4, 


99 99 




CLOSE 


#4 






OVER: 









GETNEXTWEATHER: 

'******************************* 

'GET THE NEXT WEATHER DATA BIN 
'******************************* 

NEXT j 

'******************************* 

'START CALC FOR THE NEXT MONTH 
'******************************** 

NEXT i 
/ 

CALL ENRGYCST(KWPK() , KWHTOT ( ) , BILLDEM() , BILL ( ) , YRBILL, 

PMPKWH ( ) , FANKWH(), CMPKWH ( ) , CITY) 

CALL SYSCOST ( PMPHP , FANHP, TWRCST, diameter, CITY, CAPCST) 

/ 

END 



SUB BRENT (TWIN1, TWIN2 , DIFF1, DIFF2 , TWOUT1, TWOUT2 , NTU, fda, 
MDOTW, MDOTA, TDB, TWB, DELTATCOND1 , CTR2%, TCI, TC2 , CITY) 

ITMAX =50 
TOL = .0001 
EPS = 1 * 10 A -7 
DIFF3 = DIFF2 
TWOUT3 = TWOUT2 
FOR ITER = 1 TO ITMAX 
CTR2% = CTR2% + 1 
IF (DIFF2 * DIFF3 ) > 0 THEN 
TWIN3 = TWIN1 
DIFF3 = DIFF1 
TWOUT3 = TWOUT1 
TWIN4 = TWIN2 - TWIN1 
TWINS = TWIN4 
END IF 

IF ABS (DIFF3 ) < ABS (DIFF2 ) THEN 
TWIN1 = TWIN2 
TWIN2 = TWIN3 
TWIN3 = TWIN1 
DIFF1 = DIFF2 
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TW0UT1 = TWOUT2 
DIFF2 = DIFF3 
TWOUT2 = TWOUT3 
DIFF3 = DIFF1 
TWOUT3 = TWOUT1 
END IF 

TOL1 = (2 * EPS * ABS (TWIN2 ) ) + (.5 * TOL) 

XM = .5 * (TWIN3 - TWIN 2 ) 

IF (ABS (XM) <= TOL1) OR (DIFF2 = 0) THEN 
TC2 = TWIN2 
TCI = TWOUT2 
EXIT SUB 
END IF 

IF (ABS (TWINS) >= TOL1) AND (ABS (DIFF1) > ABS (DIFF2 ) ) THEN 
S = DIFF2 / DIFF1 
IF TWIN1 = TWIN3 THEN 
P = 2 * XM * S 
Q = 1 - S 
ELSE 

Q = DIFF1 / DIFF3 
R = DIFF2 / DIFF3 

P = S * (2 * XM * Q * (Q - R) - (B - A) * (R - 1) ) 

Q = (Q - 1) * (R - 1) * (S - 1) 

END IF 

IF P > 0 THEN Q = -Q 
P = ABS (P) 

MINI = 3 * XM * Q - ABS (TOL1 * Q) 

MIN2 = ABS (TWIN5 * Q) 

IF MINI < MIN2 THEN 
MIN = MINI 
ELSE 

MIN = MIN2 
END IF 

IF (2 * P) < MIN THEN 
TWINS = TWIN4 
TWIN4 = P / Q 
ELSE 

TWIN4 = XM 
TWIN5 = TWIN4 
END IF 
ELSE 

TWIN4 = XM 
TWINS = TWIN4 
END IF 

TWIN1 = TWIN2 
DIFF1 = DIFF2 
TWOUT1 = TWOUT2 
IF (ABS (TWIN4 ) > TOL1) THEN 
TWIN2 = TWIN2 + TWIN4 
ELSE 

SELECT CASE XM 
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CASE IS > 0 
SIGN = 1 
CASE 0 

SIGN = 1 
CASE IS < 0 
SIGN = -1 
END SELECT 
T0L2 = T0L1 * SIGN 
TWIN2 = TWIN2 + TOL2 
END IF 

CALL EFFECTIVE (NTU, TWIN2, MDOTW, MDOTA, TWB, TDB , TWB, TWOUT2 ) 
DELTA2 = TWIN2 - TWOUT2 
DIFF2 = DELTATCOND1 - DELTA2 
NEXT ITER 

SELECT CASE CITY 
CASE 1 

OPEN "B: ATLEFF.DAT" FOR APPEND AS #1 
CASE 2 

OPEN "B: LAEFF.DAT" FOR APPEND AS #1 
CASE 3 

OPEN "B: NYEFF.DAT" FOR APPEND AS #1 
END SELECT 

PRINT #1, "BRENT METHOD EXCEEDED MAX ITERATIONS" 

TC2 = TWIN2 
TCI = TWOUT2 
CLOSE #1 

PRINT "BRENT METHOD EXCEEDED MAX ITERATIONS" 

PRINT "MONTH ="; i, "BIN ="; TDB 
END SUB 



SUB EFFECTIVE (NTU, TWIN, MDOTW, MDOTA, TWB, TDB, EXPTWO, TW2) 

DELTA3 =10 

Ttemp = EXPTWO 

DO WHILE ABS (DELTA3 ) > .0001 

f 

'Calculate the inlet and outlet saturation enthalpy of the 
'saturated air. 

f 

CALL Hairs (Ttemp, WSAT1, hswo) 

CALL Hairs (TWIN, WSAT1, hswi) 

i 

'Calculate the saturation specific heat and the m* parameter as 
'defined by Braun. 

f 

Cs = (hswi - hswo) / (TWIN - Ttemp) 
m = MDOTA / (MDOTW * (1 / Cs) ) 

eff = (1 - EXP (-NTU * (1 - m))) / (1 - m * EXP(-NTU * (1 - m) ) ) 

r 

'Calculate the humidity ratio of the air-vapor mixture. 
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CALL Hairs (TWB, WSAT, H) 

CALL Hairm(TWB, TDB, WSAT, 0MEG1 , hai) 
hao = hai + eff * (hswi - hai) 

Tref =32 

TW2 = ( (MDOTW * (TWIN - Tref) - MDOTA * (hao - hai)) / (MDOTW * 1)) 
+ Tref 

DELTA3 = TW2 - Ttemp 
Ttemp = TW2 
LOOP 
END SUB 



SUB ENRGYCST (KWPK(), KWHTOT ( ) , BILLDEM ( ) , BILL ( ) , YRBILL, 
PMPKWH ( ) , FANKWH ( ) , CMPKWH ( ) , CITY) 

TOTKWH = 0 

FOR 1 = 1 TO 12 

TOTKWH = TOTKWH + KWHTOT (1) 

KWHTOT ( 1) = KWHTOT ( 1 ) + 360585 
KWPK(l) = KWPK(l) + 946 
NEXT 1 

SELECT CASE CITY 
CASE 1 

GOTO ATLANTA 
CASE 2 

GOTO LOSANGELES 
CASE 3 

GOTO NEWYORK 
END SELECT 
ATLANTA: 

'************************************************************* 

' CALCULATION OF MONTHLY ENERGY COST IN ATLANTA 

' Georgia Power schedule PL_8_GS 

' General Commercial Service 

' Current as of June 1989 

/ 

KWDEM = 5 
FOR i = 6 TO 9 

KWSMDEM = .95 * KWPK(i) 

IF KWSMDEM > KWDEM THEN KWDEM = KWSMDEM 
NEXT i 

FOR i = 1 TO 12 

IF (i > 5) AND (i < 10) THEN GOTO ATLNOTWNTR 
KWWNTR = .6 * KWPK(i) 

IF KWWNTR > KWDEM THEN KWDEM = KWWNTR 
ATLNOTWNTR: 

NEXT i 

FOR i = 6 TO 9 

IF KWPK(i) > KWDEM THEN 
BILLDEM ( i) = KWPK(i) 

ELSE 
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BILLDEM(i) = KWDEM 
END IF 
NEXT i 

FOR i = 1 TO 5 

BILLDEM (i) = KWDEM 
NEXT i 

FOR i = 10 TO 12 

BILLDEM (i) = KWDEM 
NEXT i 

******************************************************** 
' CALCULATION OF THE MONTHLY ENERGY BILL FOR ATLANTA 

f******************************************************* 

YRBILL = 0 
PMPENGY = 0 
FANENGY = 0 
CMPENGY = 0 
FOR i = 1 TO 12 

BLOKHRS = KWHTOT(i) / BILLDEM (i) 

IF (BLOKHRS <= 200) THEN 
GOTO BASEI 

THEN 



ELSEIF (BLOKHRS <= 400) 
GOTO BASEI I 

ELSEIF (BLOKHRS <= 600) 

GOTO BASEIII 

ELSE 

GOTO BASEIV 
END IF 



THEN 



BASEI : 

IF KWHTOT(i) <= 3000 THEN 

COST1 = (.10185 * KWHTOT ( i ) ) +13.5 
ELSEIF KWHTOT (i) <= 10000 THEN 

COST1 = (.09282 * KWHTOT ( i ) ) + 40.59 
ELSEIF KWHTOT (i) <= 200000 THEN 

COST1 = (.07918 * KWHTOT ( i ) ) + 176.99 
ELSE 

COST1 = (.05972 * KWHTOT ( i ) ) + 4068.99 
END IF 

GOTO ATLSUMBILL 



BASEI I: 

IF (BILLDEM (i) >= 5) AND (BILLDEM (i) <= 15) THEN 

COST1 = (18.314 * BILLDEM (i) ) + (.01028 * KWHTOT ( i ) ) + 13.5 
ELSEIF (BILLDEM (i) <= 50) THEN 

COST1 = (16.508 * BILLDEM (1) ) + (.01028 * KWHTOT ( i) ) + 40.59 
ELSEIF (BILLDEM (i) <= 1000) THEN 

COST1 = (13.78 * BILLDEM (i) ) + (.01028 * KWHTOT ( i ) ) + 176.99 
ELSE 

COST1 = (9.888 * BILLDEM (i) ) + (.01028 * KWHTOT ( i ) ) + 4068.99 
END IF 

GOTO ATLSUMBILL 
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BASEIII : 

IF (BILLDEM (i) >= 5) AND (BILLDEM (i) <= 15) THEN 

C0ST1 = (19.938 * BILLDEM (i) ) + (.00622 * KWHTOT ( i ) ) +13.5 
ELSEIF BILLDEM (i) <= 50 THEN 

COST1 = (18.132 * BILLDEM (i) ) + (.00622 * KWHTOT (i) ) + 40.59 
ELSEIF BILLDEM (i) <= 1000 THEN 

COST1 = (15.404 * BILLDEM (i) ) + (.00622 * KWHTOT ( i ) ) + 176.99 
ELSE 

COST1 = (11.512 * BILLDEM (i) ) + (.00622 * KWHTOT ( i ) ) + 4068.99 
END IF 

GOTO ATLSUMBILL 

r 

BASEIV: 

IF (BILLDEM (i) >= 5) AND (BILLDEM (i) <= 15) THEN 

COST1 = (20.46 * BILLDEM (i) ) + (.00535 * KWHTOT (i) ) +13.5 
ELSEIF BILLDEM (i) <= 50 THEN 

COST1 = (18.654 * BILLDEM (i) ) + (.00535 * KWHTOT ( i) ) + 40.59 
ELSEIF BILLDEM (i) <= 1000 THEN 

COST1 = (15.926 * BILLDEM (i) ) + (.00535 * KWHTOT (i) ) + 176.99 
ELSE 

COST1 = (12.034 * BILLDEM (i) ) + (.00535 * KWHTOT ( i ) ) + 4068.99 
END IF 

' ****** CALC. OF THE TOTAL BILL ******** 

ATLSUMBILL: 

FCADJ = .016896 * KWHTOT (i) 

TAX = .05 * (COST1 + FCADJ) 

BILL ( i) = COST1 + FCADJ + TAX 
• *****CHECK FOR MIN. BILL************ 

CHKKW = BILLDEM ( i) - 30 
IF CHKKW > 0 THEN 

MINBILL = (13.5 + (7.5 * CHKKW) + FCADJ) * 1.05 
ELSE 

MINBILL = (13.5 + FCADJ) * 1.05 
END IF 

IF MINBILL >= BILL(i) THEN BILL(i) = MINBILL 
YRBILL = YRBILL + BILL(i) 

PMPENGY = PMPENGY + PMPKWH(i) 

FANENGY = FANENGY + FANKWH(i) 

CMPENGY = CMPENGY + CMPKWH(i) 

NEXT i 



FOR j = 1 TO 12 

OPEN "B: ATLEFF.DAT" FOR APPEND AS 



PRINT #4, 

PRINT #4, 

PRINT #4, 

PRINT #4, 

PRINT #4, 

PRINT #4, 

IF j = 12 THEN 



"THE MAX MONTHLY KW 
"THE BILLING DEMAND 
"THE MONTHLY KWH 
"THE MONTHLY BILL 

ii ii 



/ 

—ii 

/ 

—ii 

/ 

= •• 



#4 

KWPK ( j ) 
BILLDEM ( j ) 
KWHTOT ( j) 
BILL ( j ) 
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PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
END IF 
CLOSE #4 
NEXT j 
EXIT SUB 



"THE YEARLY TOTAL ($) YRBILL 

"THE PUMP ENERGY (KWH) =", PMPENGY 
"THE FAN ENERGY (KWH) =", FANENGY 
"THE CHILLER ENERGY (KWH) =", CMPENGY 

"THE TOTALSYSTEM ENERGY (KWH) =", TOTKWH 

•• •• 

•• •• 



LOSANGELES : 

'**************************************************************** 
'**************************************************************** 
' CALCULATION OF MONTHLY ENERGY COSTS IN LOS ANGELES 

' Souther California Edison Schedule GS-2 

' General Commercial - Large 

' Eff. Date 1 June 1988 

' Current as of June 1989 

/ 

' ** IN LA THE MONTHLY DEMAND IS STRICTLY DET. BY THE CURRENT 

' MONTH ' S PERFORMANCE. THEREFORE THE DEMAND COST IS AS FOLLOWS: 

'**************************************************************** 

YRBILL = 0 
PMPENGY = 0 
FANENGY = 0 
CMPENGY = 0 
FOR i = 1 TO 12 

BILLDEM ( i) = KWPK(i) 

DEMCOST = 0 

IF (i > 5) AND (i < 10) GOTO LASUMMER 
MIN =31.1 

DEMCOST = 2.7 * BILLDEM (i) 

IF KWHTOT(i) <= 180000 THEN 

ENGCOST = .07197 * KWHTOT(i) 

ELSEIF KWHTOT(i) > 180000 THEN 

ENGCOST = 12954.6 + ( (KWHTOT (i) - 180000) * .04973) 

ELSE 
END IF 

GOTO LASUMBILL 
LASUMMER: 

MIN =31.1 

DEMCOST = 8.6 * BILLDEM (i) 

IF KWHTOT (i) <= 180000 THEN 

ENGCOST = .07833 * KWHTOT ( i ) 

ELSEIF KWHTOT (i) > 180000 THEN 

ENGCOST = 14099.4 + ( (KWHTOT (i) - 180000) * .05247) 

ELSE 
END IF 
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LASUMBILL: 

C0ST1 = MIN + DEMCOST + ENGCOST 
SURCOST = .05 * COST1 
BILL ( i) = COST1 + SURCOST 
YRBILL = YRBILL + BILL (i) 
PMPENGY = PMPENGY + PMPKWH(i) 
FANENGY = FANENGY + FANKWH(i) 
CMP ENG Y = CMP ENG Y + CMPKWH(i) 
NEXT i 



FOR j = 1 TO 12 

OPEN "B: LAEFF.DAT" FOR APPEND AS #4 
PRINT #4, "THE MAX MONTHLY KW IS =", 

PRINT #4, "THE BILLING DEMAND =", 

"THE MONTHLY KWH =" , 

"THE MONTHLY BILL =" , 



PRINT #4, 

PRINT #4, 

PRINT #4, 

PRINT #4, 

IF j = 12 THEN 



PRINT 


#4, 


"THE 


PRINT 


#4, 


"THE 


PRINT 


#4, 


"THE 


PRINT 


#4, 


"THE 


PRINT 


#4, 


"THE 


PRINT 


#4, 


•i ii 


PRINT 


#4, 


ii n 



($) =' 



KWPK ( j ) 
BILLDEM ( j ) 
KWHTOT ( j ) 
BILL( j) 



YRBILL 
PMPENGY 
' , FANENGY 
=", CMPENGY 



END IF 
CLOSE #4 
NEXT j 
EXIT SUB 



NEWYORK : 

'**************************************************************** 
'**************************************************************** 
' CALCULATION OF MONTHLY ENERGY COSTS IN NEW YORK (CON. ED.) 

' Con Edison Rate Schedule No. 9 

' General Commercial - Large 

' Eff. Date 1 Apr. 1987 

' Current as of July 1989 

' ** IN NY THE MONTHLY DEMAND IS STRICTLY DET. BY THE CURRENT 

' MONTHS PERFORMANCE. THEREFORE THE DEMAND COST IS AS FOLLOWS: 

'**************************************************************** 

YRBILL = 0 
PMPENGY = 0 
FANENGY = 0 
CMPENGY = 0 
FOR i = 1 TO 12 

BILLDEM ( i) = KWPK(i) 

DEMCOST = 0 

IF (i > 5) AND (i < 10) GOTO NYSUMMER 
MIN = 74.18 
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IF (BILLDEM(i) > 5) 
DEMCOST = 15.86 * 
ELSEIF BILLDEM (i) > 
DEMCOST = 14274 + 
ELSE 
END IF 

GOTO NYENGYCOST 
NYSUMMER: 

MIN = 96.98 
IF (BILLDEM (i) > 5) 
DEMCOST = 20.36 * 
ELSEIF BILLDEM (i) > 
DEMCOST = 18324 + 
ELSE 
END IF 



AND (BILLDEM (i) <= 900) THEN 
BILLDEM (i) 

900 THEN 

( (BILLDEM (i) - 900) * 13.68) 



AND (BILLDEM (i) <= 900) THEN 
BILLDEM (i) 

900 THEN 

( (BILLDEM(i) - 900) * 18.18) 



************************************************************ 



' CALCULATION OF THE ENERGY CHARGES 
'*********************************************************** 



NYENGYCOST: 

IF KWHTOT(i) <= 15000 THEN 
ENGCOST = .0562 * KWHTOT(i) 

ELSEIF KWHTOT(i) > 15000 THEN 

ENGCOST = 843! + ( (KWHTOT (i) - 15000) * .0523) 
ELSE 
END IF 



COST1 = MIN + DEMCOST + ENGCOST 
MAXBILL = .3207 * KWHTOT ( i ) 

IF COST1 > MAXBILL THEN COST1 = MAXBILL 
MINBILL = MIN 

IF COST1 < MINBILL THEN COST1 = MINBILL 

/ 

FCADJ = .004414 * KWHTOT (i) 

TAX = .0825 * COST1 

BILL (i) = COST1 + FCADJ + TAX 

YRBILL = YRBILL + BILL (i) 

PMPENGY = PMPENGY + PMPKWH(i) 

FANENGY = FANENGY + FANKWH(i) 

CMPENGY = CMP ENG Y + CMPKWH(i) 

NEXT i 



FOR j = 1 TO 12 

OPEN "B: NYEFF.DAT" FOR APPEND AS #4 
PRINT #4, "THE MAX MONTHLY KW IS =", 
PRINT #4, "THE BILLING DEMAND =", 

PRINT #4, "THE MONTHLY KWH =" , 

PRINT #4, "THE MONTHLY BILL =" , 

PRINT #4, "" 

PRINT #4, "" 

IF j = 12 THEN 

PRINT #4, "THE YEARLY TOTAL ($) 



KWPK(j) 
BILLDEM (j ) 
KWHTOT ( j) 
BILL(j) 



=", YRBILL 
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PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
PRINT #4, 
END IF 
CLOSE #4 
NEXT j 



"THE PUMP ENERGY (KWH) =", PMPENGY 
"THE FAN ENERGY (KWH) =", FANENGY 
"THE CHILLER ENERGY (KWH) =", CMPENGY 
"THE TOTALSYSTEM ENERGY (KWH) =", TOTKWH 

•• ii 

•• ii 



END SUB 



SUB FLOW (GPM, CFM, TDB, MDOTW, MDOTA) 

MDOTW = GPM * 8.32124 * 60! 

MDOTA = CFM * (28.97 / (.73 * (TDB + 459.67))) * 60! 
END SUB 



SUB Hairm (T2, Tl, W2, OMEGMIXTURE, HAIRMIX) 

HWB = 1061 + .444 * T2 
HDB = 1061 + .444 * Tl 
HFLUID = T2 - 32 

OMEGMIXTURE = ((.24 * (T2 - Tl) ) + (W2 * (HWB - HFLUID))) / (HDB - 
HFLUID) 

HAIRMIX = (.24 * Tl) + OMEGMIXTURE * (1061 + .444 * Tl) 

END SUB 



SUB Hairs (TF, WS, HSAT) 

TK = (TF - 32) * (5 / 9) + 273.15 

LNPWS = -5800.2206# / TK + 1.3914993# - .04860239# * TK 

LNPWS = LNPWS + .000041764768# * (TK * 2) - .000000014452093# * (TK 

A 3) + 6.5459673# * LOG (TK) 

PWS = EXP (LNPWS) 

WS = .62198 * (PWS / (101325 - PWS)) 

HSAT = .24 * TF + WS * (1061 + .444 * TF) 

END SUB 



SUB SYSCOST ( PMPHP , FANHP, TWRCST, diameter, CITY, CAPCST) 
**************************************************************** 

' THIS SUB IS USED TO CALCULATE THE CAPITAL COST INVOLVED 
' IN THE CONDENSER SIDE EQUIPMENT INSTALLATION FOR THOSE 

' ITEMS WHICH HAVE A DIFFERENTIAL COST FOR THE SYSTEM MODS 

' DETERMINE THE COST OF THE PUMP INCLUDING ELECT INSTALLATION 
***************************************************************** 
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SELECT CASE 


PMPHP 


CASE 0 


TO 5 


PMPCST 


= 2690 


VFDCST 


= 5183 


CASE IS 


<=7.5 


PMPCST 


= 2740 


VFDCST 


= 5183 


CASE IS 


<= 10 


PMPCST 


= 3275 


VFDCST 


= 5183 


CASE IS 


<= 15 


PMPCST 


= 3765 


VFDCST 


= 6933 


CASE IS 


<= 20 


PMPCST 


= 4365 


VFDCST 


= 8618 


CASE IS 


<= 25 


PMPCST 


= 4440 


VFDCST 


= 10738 


CASE IS 


<=3 0 


PMPCST 


= 5365 


VFDCST 


= 11938 


CASE IS 


<= 50 


PMPCST 


= 6795 


VFDCST 


= 13948 


CASE IS 


> 50 


PMPCST 


= 8720 


VFDCST 


= 18050 



END SELECT 

'**************************************************************** 

' DETERMINE THE COST OF THE ELECT INSTALLATION FOR THE TOWER FAN 
'**************************************************************** 

SELECT CASE FANHP 
CASE 0 TO 5 
FANCST = 1315 
CASE IS <= 7.5 
FANCST = 1315 
CASE IS <= 10 
FANCST = 1450 
CASE IS <= 15 
FANCST = 1815 
CASE IS <= 20 
FANCST = 1965 
CASE IS <= 25 
FANCST = 1965 
CASE IS <= 30 
FANCST = 2500 
CASE IS <= 50 
FANCST = 2750 
CASE IS > 50 
FANCST = 4085 
END SELECT 
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f 'kic'k'k'k'k'k'k’k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

' DETERMINE THE CAPITAL COST OF THE PIPING, VALVING, FITTINGS 
r *********************************************************** 

SELECT CASE diameter 
CASE 5.05 
PIPECST = 19850 
CASE 6.07 
PIPECST = 23495 
CASE 8.07 
PIPECST = 33548 
CASE 10.02 
PIPECST = 52217 
END SELECT 

• ********************************************************** 

' DETERMINE THE TOTAL CAPITAL COST 
' ********************************************************** 

TOTCST = TWRCST + PMPCST + FANCST + PIPECST + VFDCST 

OVRHD = .15 * TOTCST 

SUB1 = TOTCST + OVRHD 

PROF = .1 * SUB1 

CAPCST = SUB1 + PROF 

> ********************************************************** 

' INCLUDE THE 1989 MEANS CITY COST INDEXES 
f ********************************************************** 

SELECT CASE CITY 
CASE 1 

CTYINDX = .917 
CASE 2 

CTYINDX = 1.117 
CASE 3 

CTYINDX = 1.235 
END SELECT 

CAPCST = CAPCST * CTYINDX 
SELECT CASE CITY 
CASE 1 

OPEN "B: ATLEFF.DAT" FOR APPEND AS #4 
CASE 2 



OPEN 


"B: 


LAEFF. 


, DAT " FOR APPEND 


AS #4 


CASE 


3 










OPEN 


"B: 


NYEFF. 


.DAT" FOR APPEND 


AS #4 


END SELECT 








PRINT 


#4, 


"THE 


TOWER COST IS 


tl 

/ 


TWRCST 


PRINT 


#4, 


"THE 


PUMP COST IS 


II 

/ 


PMPCST 


PRINT 


#4, 


"THE 


FAN ELECT COST", 


FANCST 


PRINT 


#4, 


"THE 


PIPE SYS COST 


II 

/ 


PIPECST 


PRINT 


#4, 


"THE 


SUBTOTAL IS 


II 

/ 


TOTCST 


PRINT 


#4, 


"THE 


OVERHEAD 


II 

/ 


OVRHD 


PRINT 


#4, 


"THE 


PROFIT IS 


II 

/ 


PROF 


PRINT 


#4, 


"THE 


CITY INDEX IS 


II 

/ 


CTYINDX 


PRINT 


#4, 


"THE 


TOTAL COST IS 


II 

/ 


CAPCST 


CLOSE 


#4 











END SUB 
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APPENDIX B 



DATA GENERATED FROM WEBER MODEL COMPARISON 
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OUTPUT OF DELTA TEMPERATURE FOR EFFECTIVENESS AND JOYCE METHOD 



DELTA EFFECTIVENESS 


DELTA WEBER 


DELTA WEBER 


DELTA WEBER 


(DEG F) 


(DEG F) 


PLUS 15% 


MINUS 15% 


1.6 


1.4 


1.61 


1.19 


6.4 


5.9 


6.79 


5.02 


6.5 


5.9 


6.79 


5.02 


6.5 


5.9 


6.79 


5.02 


7.8 


7.3 


8.39 


6.2 


9.8 


8.9 


10.235 


7.57 


11 


10 


11.5 


8.5 


11 


10 


11.5 


8.5 


11.1 


10.4 


11.96 


8.84 


12.2 


11.2 


12.88 


9.52 


14.5 


13 


14.95 


11.05 


14 . 6 


13 


14.95 


11.05 


14 . 6 


13 


14.95 


11.05 


19.9 


17.8 


20.47 


15.13 


20 


18.1 


20.82 


15.39 


20.1 


17.8 


20.47 


15.13 


20.2 


17.8 


20.47 


15.13 


26.9 


23.4 


26.91 


19.89 


26.9 


23.4 


26.91 


19.89 


26.9 


23.4 


26.91 


19.89 


27 


23.6 


27.14 


20.06 


28.6 


25 


28.75 


21.25 


28.8 


25 


28.75 


21.25 


29 


25 


28.75 


21.25 
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APPENDIX C 



DATA GENERATED FROM JOYCE MODEL COMPARISON 
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AT750.DAT PUMP 8IZE 750 GPM JOYCE METHOD 



MONTH= 12 
START TIME 
ctrl= 7 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 



11:12:25 END TIME = 11:53:09 RUNTIME = 40:44 
Ctr2 = 92 Ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

100 

CWR(DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



TOWER 

3 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.60413 .1174005 

pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
6.151319 



NfISENT 

.5164742 

KWHBIN 

2526.094 



FANKW 

4.87361 



11 



NFTON 

.6406521 

KWHTOT 

6573.959 



TOTALKW 

101.7707 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



334021.2 

17345.86 

46195.46 

337459.7 



401001 
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AT100.DAT PUMP SIZE 1000 6PM JOYCE METHOD 



MONTH= 12 BIN =57 

START TIME = 11:55:00 END TIME = 12:37:03 RUNTIME = 42:03 

ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2481.199 



FANKW 

4.87361 



TOWER 

3 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.60413 .1174005 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6469.957 



TOTALKW 

99.96196 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



333374 

9823.913 

46195.46 

337459.7 



393479.1 
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AT125.DAT PUMP SIZE 1250 GPM JOYCE METHOD 



MONTH= 12 BIN = 57 

START TIME = 11:36:25 END TIME = 12:19:33 RUNTIME = 43:08 

ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

60 

CWR(DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.541853 



NfISENT 

.5164742 

KWHBIN 

2486.145 



FANKW 

4.87361 



TOWER 

3 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.60413 .1174005 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND.TEMP 

.6406521 84.4407 

KWHTOT 
6481.415 



TOTALKW 

100.1612 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



333380.8 

10126.37 

46195.46 

337459.7 



393781.5 



44 



AT150.DAT PUMP SIZE 1500 6PM JOYCE METHOD 



MONTH= 


12 


BIN = 57 








START TIME = 


12:25:32 


END TIME = 13: 


09:23 RUNTIME 


= 43:51 


ctrl= 7 




Ctr2 = 92 


Ctr3 = 3 






1= 12 




J= 10 








TDB (DEG 


F) 


TWB (DEG 


F) HOURS 


TOWER 




57 




52 


24.82143 


3 




GPM 




%GPM 


CFM 


TONS 




750 




50 


120000 


100 




CWS (DEG 


F) 


CWR(DEG F) 


DELTATCOND 


DELTATTOWER 


TOWERDUTY 


80 




84.03179 


4.031794 


19.60413 


. 1174005 



PIPE D (FT) Pipevelocity (FT/SEC) REd pipe frict 

.6666667 4.787051 369371.3 1.592428E-02 



pipeL total PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
5.159607 



NfISENT 

.5164742 

KWHBIN 

2501.478 



FANKW 

4.87361 



11 



NFTON 

.6406521 

KWHTOT 

6516.936 



TOTALKW 

100.779 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



333485.1 

11436.89 

46195.46 

337459.7 



395092 
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ATEF750.DAT PUMP SIZE 750 6PM EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 16:44:50 TIME END = 16:51:52 RUNTIME =7:02 
ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 490044E-02 
PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

100 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

3 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



21.74524 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PH tower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
6.151319 



NfISENT 

.5164742 

KWHBIN 

2502.91 



FANKW 

3.939574 



11 



NFTON 

.6406521 

KWHTOT 

6502.745 



TOTALKW 

100.8366 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



333296 

17345.86 

40526.96 

337296.7 



395169.5 
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ATEF100.DAT POMP SIZE 1000 6PM EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 7 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 490044E-02 
PIPE D (FT) 
.6666667 



BIN = 57 
17:22:58 



TIME END = 17:29:55 RUNTIME =6:57 



Ctr2 = 92 

J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 

84.03179 



ctr3 = 3 

HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

3 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



21.74524 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe 

PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2458.015 



FANKW 

3.939574 



PHtower 



11 



21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6398.744 



TOTALKW 

99.02793 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



332648.8 

9823.913 

40526.96 

337296.7 



387647.6 
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ATEF125.DAT PUMP 8IZE 1250 6PM EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 7 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 490044E-02 
PIPE D (FT) 

. 6666667 



BIN = 57 



17:32:44 TIME END = 17:40:57 RUNTIME =8:13 
Ctr2 = 92 Ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

60 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

3 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



21.74524 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1. 592428E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
. 9074562 
HOURSSUM 
57.5 
PUMPKW 
4.541853 



NfISENT 

.5164742 

KWHBIN 

2462.961 



FANKW 

3.939574 



11 



NFTON 
. 6406521 
KWHTOT 
6410.201 



TOTALKW 

99.22719 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) 

THE CHILLER ENERGY (KWH) = 
THE TOTALS Y STEM ENERGY (KWH) 



332655.6 
10126.37 
40526.96 

337296.7 



387950 
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ATEF150.DAT PUMP 8IZE 1500 6PM EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 17:42:20 TIME END = 17:51:43 RUNTIME = 9:23 
ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 490044E-02 
PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 


HOURS 


TOWER 


52 


24.82143 


3 


%GPM 


CFM 


TONS 


50 


120000 


100 


CWR(DEG F) 


DELTATCOND 


DELTATTOWER TOWERDUTY 


84.03179 


4.031794 21.74524 


Pipevelocity (FT/SEC) REd 


pipe frict 


4.787051 


369371.3 


1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

5.159607 



NfISENT 

.5164742 

KWHBIN 

2478.294 



FANKW 

3.939574 



11 



NFTON 

.6406521 

KWHTOT 

6445.722 



TOTALKW 

99.84494 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



332759.9 

11436.89 

40526.96 

337296.7 



389260.5 



49 



ATT0W2.DAT TOWER # 2 JOYCE METHOD 



MONTH= 12 
START TIME 
ctrl= 7 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 


10:02:55 


END TIME = 


10:43:28 RUNTIME = 40:33 


Ctr2 = 90 


Ctr3 = 3 




J= 10 


TWB ( DEG F) 


HOURS 


TOWER 


52 


24.82143 


2 


%GPM 


CFM 


TONS 


75 


120000 


100 


CWR(DEG F) 


DELTATCOND 


DELTATTOWER TOWERDUTY 


84.03179 


4.031794 


19.00449 .1246106 


Pipevelocity (FT/SEC) REd 


pipe frict 


4.787051 


369371.3 


1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074574 

HOURS 

24.82143 

CHILLERKW 

90.74574 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2541.575 



FANKW 

7.306037 



11 



NFTON 

.6406521 

KWHTOT 

6686.729 



TOTALKW 

102.3944 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



336101.8 
9823.913 
67227.36 

337547.8 



414599.1 



50 



AT TOW 5 • DAT TOWER # 5 JOYCE METHOD 



MONTH= 12 BIN = 57 

START TIME = 10:44:54 END TIME = 11:26:56 RUNTIME = 42:02 

ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 

REd 

369371.3 



Pipevelocity (FT/SEC) 
4.787051 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2481.199 



FANKW 

4.87361 



TOWER 

5 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.60413 .1174005 

pipe frict 
1.592428E-02 
PHtower 



11 



21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6469.957 



TOTALKW 

99.96196 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



333374 

9823.913 

46195.46 

337459.7 



393479.1 



51 



ATT0W9.DAT TOWER § 9 JOYCE METHOD 



MONTH= 12 
START TIME 
ctrl= 7 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 



11:29:52 END TIME = 12:13:13 RUNTIME = 43:21 
Ctr2 = 91 Ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



TOWER 

9 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.93193 .1136424 

pipe frict 
1.592428E-02 



pipeL total PHpipe 

PHtotal ( FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2405.123 



FANKW 

1.808699 



PHtower 



11 21.77203 



NFTON COND.TEMP 

.6406521 84.4407 

KWHTOT 
6187.521 



TOTALKW 

96.89704 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329604.7 

9823.913 

17438.78 

337425.3 



364687.9 



52 



ATT0W12.DAT TOWER § 12 JOYCE METHOD 



MONTH= 12 
START TIME 
ctn= 7 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 

12:21:59 END TIME = 13:02:36 RUNTIME = 40:37 



Ctr2 = 91 

J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 
84.03179 



ctr3 = 3 

HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



TOWER 

12 

TONS 

100 

DELTATTOWER TOWERDUTY 

19.7356 .1158782 

pipe frict 
1 . 592428E-02 



pipeL total PHpipe 
PHtotal(FT Water) 



PHcond 



666.6667 

PASSES 



5.666439 5.105591 

Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2397.792 



FANKW 

1.513313 



11 



NFTON 

.6406521 

KWHTOT 

6159.366 



TOTALKW 

96.60166 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329208.7 

9823.913 

14440.33 

337443.6 



361707.9 



53 



ATTOW 16. DAT TOWER # 16 JOYCE METHOD 



MONTH= 12 BIN =57 

START TIME = 13:04:09 END TIME = 13:46:03 RUNTIME = 41:54 

ctrl= 7 ctr2 = 90 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
. 6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR (DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 

REd 

369371.3 



Pipevelocity (FT/SEC) 
4.787051 



pipeL total PHpipe 

PHtotal(FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2385.629 



FANKW 

1.023324 



TOWER 

16 

TONS 

100 

DELTATTOWER TOWERDUTY 

20.19933 .1106671 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND. TEMP 

.6406521 84.4407 

KWHTOT 
6114 . 878 



TOTALKW 

96.11167 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS Y STEM ENERGY (KWH) 



328628.1 
9823.913 
10003.45 

337395.2 



357222.5 



54 



ATT0W19.DAT TOWER # 19 JOYCE METHOD 



MONTH= 12 
START TIME 
ctrl= 7 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 



14:35:43 END TIME = 15:19:20 RUNTIME = 43:37 
ctr2 = 90 ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



TOWER 

19 

TONS 

100 

DELTATTOWER TOWERDUTY 

20.77497 .104522 

pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PH total (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074577 

HOURS 

24.82143 

CHILLERKW 

90.74577 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2385.587 



FANKW 

1.021625 



11 



NFTON 

.6406521 

KWHTOT 

6115.9 



TOTALKW 

96.10999 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



328658.3 

9823.913 

10252.19 

337320.8 



357396.9 



55 



ATEF2.DAT TOWER § 2 EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 15:37:19 TIME END = 15:44:00 RUNTIME =6:41 
ctrl= 7 ctr2 = 91 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2509.533 



FANKW 

6.015156 



TOWER 

2 

TONS 

100 

DELTATTOWER TOWERDUTY 

20.96244 .1025935 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6588.154 



TOTALKW 

101.1035 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



335101.8 

9823.913 

59464.87 

337298.5 



406587.3 



56 



ATEF5.DAT TOWER # 5 EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 15:46:17 TIME END = 15:54:27 RUNTIME =8:10 
ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 490044E-02 
PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

5 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



21.74524 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 



5.666439 5.105591 

Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



2.900009 

LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 



NfISENT 

.5164742 

KWHBIN 

2458.015 



FANKW 

3.939574 



11 



NFTON 

.6406521 

KWHTOT 

6398.744 



TOTALKW 

99.02793 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



332648.8 

9823.913 

40526.96 

337296.7 



387647.6 



57 



ATEF9.DAT TOWER # 9 EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 7 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 089065E-02 
PIPE D (FT) 
.6666667 



BIN = 57 



15:57:15 TIME END = 16:06:35 RUNTIME = 9:20 
ctr2 = 92 ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR (DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

9 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



22.17689 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe 
PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074577 

HOURS 

24.82143 

CHILLERKW 

90.74577 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2396.136 



FANKW 

1.446589 



PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6159.948 



TOTALKW 

96.53496 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329324.2 

9823.913 

15236.21 

337296.9 



362357 



58 



ATEF12.DAT TOWER # 12 EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 7 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

9 . 327599E-02 
PIPE D (FT) 
.6666667 



BIN = 57 



16:28:21 TIME END = 16:35:14 RUNTIME =6:53 
ctr2 = 93 ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

12 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



21.91807 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1. 592428E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 
893100.2 
KWTON GUESS 
.9074577 
HOURS 
24.82143 
CHILLERKW 
90.74577 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2390.465 



FANKW 

1.218139 



11 



NFTON 

.6406521 

KWHTOT 

6136.881 



TOTALKW 

96.30651 



21.77203 

COND.TEMP 

84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



328978.1 

9823.913 

12646.48 

337296.8 



359767.2 



59 



ATEF16.DAT TOWER # 16 EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 7 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

8 . 771837E-02 
PIPE D (FT) 

. 6666667 



BIN = 57 
16:36:34 



TIME END = 16:44:47 



Ctr2 = 91 

J= 10 

TWB (DEG F) 
52 

%GPM 
75 

CWR (DEG F) 

84.03179 



Ctr3 = 3 

HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



RUNTIME = 8:13 



TOWER 

16 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



22.53072 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 
893100.2 
KWTONGUES S 
.9074577 
HOURS 
24.82143 
CHILLERKW 
90.74577 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2380.363 



FANKW 

.8111207 



11 



NFTON 

.6406521 

KWHTOT 

6098.738 



TOTALKW 

95.89949 



21.77203 

COND. TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



328464.3 

9823.913 

8709.041 

337297.1 



355830.1 



60 



ATEF19.DAT TOWER # 19 EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 8 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

8 . 117554E-02 
PIPE D (FT) 
.6666667 



BIN = 57 



16:46:15 TIME END = 16:55:43 RUNTIME =9:28 
Ctr2 = 104 Ctr3 = 3 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 



TOWER 

19 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



23.29736 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074578 

HOURS 

24.82143 

CHILLERKW 

90.74578 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2379.924 



FANKW 

.7934307 



11 



NFTON 

.6406521 

KWHTOT 

6098.561 



TOTALKW 

95.88181 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



328486.1 

9823.913 

8849.848 

337296.8 



355970.5 



61 



ATFM08.DAT CFM 80,000 



JOYCE METHOD 



MONTH= 12 BIN =57 

START TIME = 14:10:57 END TIME = 14:51:45 RUNTIME = 40:48 
ctrl= 6 ctr2 = 80 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 
pipeL total 



HOURS 

24.82143 

CFM 

80000 

DELTATCOND 

4.031793 



J= 10 

TWB ( DEG F) 

52 

%GPM 
75 

CWR(DEG F) 

84.03179 

Pipevelocity (FT/SEC) REd 
4.787051 369371.3 

PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

Condenservelocity 
2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074569 

HOURS 

24.82143 

CHILLERKW 

90.74569 



NfISENT 

.5164742 

KWHBIN 

2386.287 



FANKW 

1.049877 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

14.52383 .1973316 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6108.143 



TOTALKW 

96.13816 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329347.2 
9823.913 
7669.604 

341085.2 



358578.7 
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ATFM10.DAT 



CFH 100,000 



JOYCE METHOD 



MONTH= 12 BIN =57 

START TIME = 15:34:44 END TIME = 16:17:07 RUNTIME = 42:23 
ctrl= 6 ctr2 = 80 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR (DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

100000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 

PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074571 

HOURS 

24.82143 

CHILLERKW 

90.74571 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
. 9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2401.928 



FANKW 

1.680003 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

16.90012 .1539621 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6167.014 



TOTALKW 

96.76831 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329446 

9823.913 

13997.07 

338508.3 



362329.2 
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ATFM12 . DAT 



CFM 120,000 JOYCE METHOD 



MONTH= 12 BIN = 57 

START TIME = 09:38:23 END TIME = 10:19:09 RUNTIME = 40:46 
ctrl= 7 ctr2 = 91 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
. 6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 

REd 

369371.3 



Pipevelocity (FT/SEC) 
4.787051 



pipeL total PHpipe 

PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074574 

HOURS 

24.82143 

CHILLERKW 

90.74574 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2421.411 



FANKW 

2.464895 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

18.92442 .1256079 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6244.041 



TOTALKW 

97.55323 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



330303.3 
9823.913 
22643.94 

337620.4 



370088.2 
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ATFM14.DAT 



CFM 140/ 000 JOYCE METHOD 



MONTH= 12 BIN = 57 

START TIME = 09:20:36 END TIME = 10:03:05 RUNTIME = 42:29 
ctrl= 7 ctr2 = 90 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 
pipeL total 



HOURS 

24.82143 

CFM 

140000 

DELTATCOND 

4.031794 



J= 10 

TWB (DEG F) 

52 

%GPM 
75 

CWR ( DEG F) 

84.03179 

Pipevelocity (FT/SEC) REd 
4.787051 369371.3 

PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

Condenservelocity 
2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 



PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074576 

HOURS 

24.82143 

CHILLERKW 

90.74576 



NfISENT 

.5164742 

KWHBIN 

2444.862 



FANKW 

3.409683 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

20.65385 .1057865 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6341.728 



TOTALKW 

98.49804 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



331787.7 

9823.913 

34027.18 

337324.1 



381175.2 
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ATFM16.DAT CFM 160,000 


JOYCE METHOD 




MONTH = 12 


BIN = 57 








START TIME = 


10:04:59 END 


TIME = 10:48 


: 51 RUNTIME = 


43:52 


ctrl= 8 


Ctr2 = 100 


ctr3 = 3 






1= 12 


J= 10 








TDB (DEG F) 


TWB (DEG F) 


HOURS 


TOWER 




57 


52 


24.82143 


8 




GPM 


%GPM 


CFM 


TONS 




750 


75 


160000 


100 




CWS (DEG F) 


CWR ( DEG F) 


DELTATCOND 


DELTATTOWER 


TOWERDUTY 


80 


84.03179 


4.031794 22.13515 


9 . 127153E-02 










PIPE D (FT) 


Pipevelocity (FT/SEC) REd 


pipe frict 


.6666667 


4.787051 


369371.3 


1 . 592428E-02 




pipeL total 


PHpipe 


PHcond 


PHtower 




PHtotal(FT Water) 








666.6667 


5.666439 


5.105591 


11 


21.77203 


PASSES 


Condenservelocity 






2 


2.900009 








UACOND 


LMTDCOND 


NfISENT 


NFTON 


COND . TEMP 


893100.2 


1.690421 


.5164742 


.6406521 


84.4407 


KWTONGUESS 


KWTONMODEL 


KWHBIN 


KWHTOT 




. 9074579 


.9074562 


2472.463 


6462.413 




HOURS 


HOURSSUM 








24.82143 


57.5 








CHILLERKW 


PUMPKW 


FANKW 


TOTALKW 




90.74579 


4.342596 


4.521638 


99.61002 





THE 


YEARLY TOTAL ($) = 


333746.1 




THE 


PUMP ENERGY (KWH) = 


9823.913 




THE 


FAN ENERGY (KWH) = 


48544.85 




THE 


CHILLER ENERGY (KWH) = 


337296.6 




THE 


TOTALSYSTEM ENERGY (KWH) 


= 


395665.4 
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ATFM18.DAT CFM 180,000 JOYCE METHOD 



MONTH= 12 
START TIME 
Ctrl= 8 
1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



BIN = 57 



13:28:34 END TIME = 14:09:48 RUNTIME = 41:14 



Ctr2 = 94 

J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 
84.03179 



Ctr3 = 3 

HOURS 

24.82143 

CFM 

180000 

DELTATCOND 

4.031794 



Pipevelocity ( FT/SEC) REd 
4.787051 369371.3 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

23.40676 .0802768 

pipe frict 
1 . 592428E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 

893100.2 

KWTONGUESS 

.9074581 

HOURS 

24.82143 

CHILLERKW 

90.7458 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2504.454 



FANKW 

5.810467 



11 



NFTON 

.6406521 

KWHTOT 

6608.818 



TOTALKW 

100.8989 



21.77203 

COND . TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



336225.5 

9823.913 

66800.08 

337297.3 



413921.3 
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APMEF08.DAT CFM 80,000 



EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 14:16:25 TIME END = 14:24:02 RUNTIME =7:37 
ctrl= 6 ctr2 = 85 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
. 6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

80000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 

PHtotal (FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.907457 

HOURS 

24.82143 

CHILLERKW 

90.7457 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 

1.690421 

KWTONMODEL 

.9074562 

HOURSSUM 

57.5 

PUMPKW 

4.342596 



NfISENT 

.5164742 

KWHBIN 

2383.514 



FANKW 

.9381654 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

15.58475 .1763346 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6099.48 



TOTALKW 

96.02646 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS Y STEM ENERGY (KWH) 



328801.4 
9823.913 
7206.588 

339184.5 



356215 
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AFMEF10.DAT CFM 100,000 



EFFECTIVENESS METHOD 



MONTH= 12 BIN =57 

TIME START = 14:28:05 TIME END = 14:36:57 RUNTIME =8:52 
ctrl= 6 ctr2 = 81 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

100000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe 

PHtotal(FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074572 

HOURS 

24.82143 

CHILLERKW 

90.74572 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2396.03 



FANKW 

1.442381 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

18.4128 .1321855 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6148.765 



TOTALKW 

96.53069 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALS YSTEM ENERGY (KWH) 



329002.8 

9823.913 

12707.91 

337642.3 



360174.1 
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AFMEF12.DAT 



CFM 120, 000 



EFFECTIVENESS METHOD 



MONTH = 12 BIN = 57 

TIME START = 12:46:37 TIME END = 12:53:25 RUNTIME =6:48 
ctrl= 7 ctr2 = 92 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 
pipeL total 



HOURS 

24.82143 

CFM 

120000 

DELTATCOND 

4.031794 



J= 10 

TWB (DEG F) 

52 

%GPM 
75 

CWR(DEG F) 

84.03179 

Pipevelocity (FT/SEC) REd 
4.787051 369371.3 

PHpipe PHcond 

PHtotal(FT Water) 

666.6667 5.666439 5.105591 

Condenservelocity 
2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURSSUM 
57.5 
PUMPKW 
4.342596 



PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074575 

HOURS 

24.82143 

CHILLERKW 

90.74575 



NfISENT 

.5164742 

KWHBIN 

2410.719 



FANKW 

2.034147 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

20.85859 .1036575 

pipe frict 
1 . 592428E-02 
PHtower 



11 



21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6211.164 



TOTALKW 

97.12249 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



329947.8 

9823.913 

20089.43 

337312.3 



367225.7 



70 



AFMEF14.DAT CFH 140,000 EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 12:54:56 TIME END = 13:03:20 RUNTIME =8:24 
ctrl= 8 ctr2 = 104 ctr3 = 3 



1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

8 . 405238E-02 
PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR(DEG F) 

84.03179 



HOURS 

24.82143 

CFM 

140000 

DELTATCOND 



TOWER 

8 

TONS 

100 

DELTATTOWER 



TOWERDUTY 



4.031794 



22.95394 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipe frict 
1 . 5924 28E-02 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

PASSES Condenservelocity 



PHtower 



2 

UACOND 
893100.2 
KWTONGUES S 
.9074579 
HOURS 
24.82143 
CHILLERKW 
90.74579 



2.900009 
LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2427.475 



FANKW 

2.709155 



11 



NFTON 

.6406521 

KWHTOT 

6288.67 



TOTALKW 

97.79754 



21.77203 

COND. TEMP 
84.4407 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



331261.7 

9823.913 

29775.94 

337297 



376896.8 
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AFMEF16.DAT CFM 160,000 



EFFECTIVENESS METHOD 



MONTH= 12 BIN = 57 

TIME START = 11:06:16 TIME END = 11:13:18 RUNTIME =7:02 
ctrl= 8 ctr2 = 105 ctr3 = 3 



1 = 12 

TDB (DEG F) 
57 
GPM 
750 

CWS (DEG F) 
80 

PIPE D (FT) 
.6666667 



J= 10 

TWB (DEG F) 
52 

%GPM 

75 

CWR ( DEG F) 
84.03179 



HOURS 

24.82143 

CFM 

160000 

DELTATCOND 

4.031794 



Pipevelocity (FT/SEC) REd 
4.787051 369371.3 



pipeL total PHpipe PHcond 

PHtotal (FT Water) 

666.6667 5.666439 5.105591 

Condenservelocity 
2.900009 

LMTDCOND NfISENT 

1.690421 .5164742 

KWTONMODEL KWHBIN 

.9074562 2446.345 

HOURS SUM 
57.5 

PUMPKW FANKW 

4.342596 3.469392 



PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074581 

HOURS 

24.82143 

CHILLERKW 

90.7458 



TOWER 

8 

TONS 

100 

DELTATTOWER TOWERDUTY 

24.73064 .0700314 

pipe frict 
1 . 592428E-02 
PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6383.272 



TOTALKW 

98.55779 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



332982.5 

9823.913 

42322.18 

337297.3 



389443.3 
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AFMEF18 . DAT 



CFM 180,000 



EFFECTIVENESS METHOD 



MONTH= 12 
TIME START = 
ctrl= 9 
1 = 12 

TDB (DEG F) 

57 

GPM 

750 

CWS (DEG F) 

80 

5 . 973898E-02 
PIPE D (FT) 
.6666667 



BIN = 57 

12:18:54 TIME END = 12:27:39 RUNTIME =8:45 



ctr2 = 116 


ctr3 = 3 




J= 10 

TWB (DEG F) 


HOURS 


TOWER 


52 


24.82143 


8 


%GPM 


CFM 


TONS 


75 


180000 


100 


CWR(DEG F) 


DELTATCOND 


DELTATTOWER TOWERDUTY 


84.03179 


4.031795 26.22048 


Pipevelocity (FT/SEC) REd 


pipe frict 


4.787051 


369371.3 


1 . 592428E-02 



pipeL total PHpipe 
PHtotal(FT Water) 



PHcond 



666.6667 

PASSES 

2 

UACOND 

893100.2 

KWTONGUESS 

.9074584 

HOURS 

24.82143 

CHILLERKW 

90.74583 



5.666439 5.105591 

Condenservelocity 
2.900009 



LMTDCOND 
1.690421 
KWTONMODEL 
.9074562 
HOURS SUM 
57.5 
PUMPKW 
4.342596 



NfISENT 

.5164742 

KWHBIN 

2467.557 



FANKW 

4.323931 



PHtower 



11 21.77203 



NFTON COND . TEMP 

.6406521 84.4407 

KWHTOT 
6497.896 



TOTALKW 

99.41236 



THE YEARLY TOTAL ($) = 

THE PUMP ENERGY (KWH) = 

THE FAN ENERGY (KWH) = 

THE CHILLER ENERGY (KWH) = 
THE TOTALSYSTEM ENERGY (KWH) 



335188.1 

9823.913 

58278.78 

337297.4 



405400.2 
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Iota 





r)-v 



Thesis 
R3755 
c . 1 



Thesis 
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Rice 

Employing an effective- 
ness model for calcula- 
ting cooling tower per- 
formance. 



